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Effects of collisions on conservation laws for toroidal plasmas are investigated based on the gyrokinetic field 
theory. Associating the collisional system with a corresponding collisionless system at a given time such that 
the two systems have the same distribution functions and electromagnetic fields instantaneously, it is shown 
how the collisionless conservation laws derived from Noether’s theorem are modified by the collision term. 
Effects of the external source term added into the gyrokinetic equation can be formulated similarly with 
the collisional effects. Particle, energy, and toroidal momentum balance equations including collisional and 
turbulent transport fluxes are systematically derived using a novel gyrokinetic collision operator, by which the 
collisional change rates of energy and canonical toroidal angular momentum per unit volume in the gyrocenter 
space can be given in the conservative forms. The ensemble-averaged transport equations of particles, energy, 
and toroidal momentum given in the present work are shown to include classical, neoclassical, and turbulent 
transport fluxes which agree with those derived from conventional recursive formulations. 
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I. INTRODUCTION 

Gyrokinetic theories and simulations are powerful 
means to investigate microinstabilities and turbulent 
transport processes in magnetically confined plasmasii^— 
Originally, gyrokinetic equations are derived by recur¬ 
sive techniques combined with the WKB or ballooning 
representationj^^— On the other hand, modern deriva¬ 
tions of the gyrokinetic equations are based on the La- 
grangian and/or Hamiltonian formulations^ii in which 
conservation laws for the phase-space volume and the 
magnetic moment are automatically ensured by Liou- 
ville’s theorem and Noether’s theorem, respectively^^ 
Besides, conservation of the total energy and momen¬ 
tum is naturally obtained in the gyrokinetic field the¬ 
ory, where all governing equations for the distribution 
functions and the electromagnetic fields are derived from 
the Lagrangian which describes the whole system consist¬ 
ing of particles and fieldsji^i^— A subtle point regarding 
the Lagrangian/Hamiltonian gyrokinetic formulations is 
that they basically treat collisionless systems so that 
Noether’s theorem and conservation laws do not hold di¬ 
rectly for collisional systems. In this paper, we examine 
how the collision and external source terms added into 
the gyrokinetic equations influence the conservation laws 
derived from Noether’s theorem in the gyrokinetic field 
theory for collisionless systems. 

For a given collisional kinetic system, we can imag¬ 
ine a corresponding collisionless kinetic system such that 
the two systems have the same distribution functions 
and electromagnetic fields instantaneously. As an exam¬ 
ple of two such systems, the Boltzmann-Poisson-Ampere 
system and the Vlasov-Poisson-Ampere system are con¬ 
sidered in Sec. II, where we express the variation of 
the action integral for the latter collisionless system in 
terms of the distribution functions and the electromag¬ 


netic fields for the former collisional system to show 
how the conservation laws derived from Noether’s theo¬ 
rem in the collisionless system are modified in the col¬ 
lisional system with external sources of particles, en¬ 
ergy, and momentum. There, we confirm the natural 
result that, when adding no external sources but only 
the collision term that conserves the energy and momen¬ 
tum, the energy and momentum conservation laws for 
the Boltzmann-Poisson-Ampere system take the same 
forms as those for the Vlasov-Poisson-Ampere system. 
The above-mentioned procedures are repeated in Sec. Ill 
to treat the collisional and collisionless gyrokinetic sys¬ 
tems. In our previous work^ using the gyrokinetic 
Vlasov-Poisson-Ampere system of equations, conserva¬ 
tion laws of particles, energy, and toroidal angular mo¬ 
mentum are obtained for collisionless toroidal plasmas, 
in which the slow temporal variation of the background 
magnetic field is taken into account in order to enable 
self-consistent treatment of physical processes on trans¬ 
port time scales. Based on these results, the particle, en¬ 
ergy, and toroidal angular momentum balance equations 
for the collisional plasma are derived from the gyroki¬ 
netic Boltzmann-Poisson-Ampere system of equations in 
Secs. IV and V. In Sec. VI, it is shown by taking the en¬ 
semble average of these balance equations that the par¬ 
ticle, energy, and toroidal angular momentum transport 
fluxes are given by the sum of the conventional expres¬ 
sions of the classical, neoclassical, and turbulent trans¬ 
port fluxes to the lowest order in the normalized gyro- 
radius parameter. Conclusions are given in Sec. VII and 
formulas for transformation from particle to gyrocenter 
coordinates are presented in Appendix A. 

Regarding the collision operator for the gyrokinetic 
equation, several works have been done, which take ac¬ 
count of finite-gyroradius effects to modify the Landau 
collision operator defined in the particle coordinates 
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The relation of the collision operator in the gyrocenter co¬ 
ordinates to that in the particle coordinates is explained 
in Appendix B. The Landau operator for Coulomb col¬ 
lisions conserves particles’ number, kinetic energy, and 
momentum locally at the particle position although, in 
the gyrocenter position space, collisions induce transport 
fluxes of particles, energy, and momentum. Besides, it is 
emphasized in this work that the collisional change rates 
of the gyrocenter Hamiltonian (which includes not only 
the kinetic energy but also the potential energy) and of 
the canonical momentum (instead of the kinetic momen¬ 
tum) per unit volume in the gyrocenter space should take 
the conservative (or divergence) forms in order to prop¬ 
erly derive the energy and momentum conservation laws 
for the collisional gyrokinetic system. The approximate 
collision operator which keeps these conservation proper¬ 
ties of the gyrocenter energy and canonical toroidal an¬ 
gular momentum is shown in Appendix C. [It is noted 
that another form of the gyrokinetic collision operator, 
which satisfies the energy and momentum conservation 
laws, has recently been presented by Burby et al^] Ap¬ 
pendix D is given to describe how to derive the formula 
for the toroidal angular momentum transport flux due to 
the collision term. 


II. BOLTZMANN-POISSON-AMPERE SYSTEM 


In this section, conservation laws are investigated 
for the Boltzmann-Poisson-Ampere system of equations 
which provide the basis of approximate description by the 
collisional electromagnetic gyrokinetic system of equa¬ 
tions for strongly magnetized plasmas considered in 
Secs. Ill-VI. Time evolution of the distribution function 
/o(x, V, t) for particle species a is described by the Boltz¬ 
mann kinetic equation. 


+ vV + -^ 
rUa 

= A:a(x, V,t), 


E(x, f) -I—V X B(x,t) 
c 


A 

dv 


( 1 ) 


where ICa{x.,'v,t) denotes the rate of change in the dis¬ 
tribution function fa due to Coulomb collisions and it 
may also include other parts representing external par¬ 
ticle, momentum, and/or energy sources if any. The 
electromagnetic fields E(x, t) and B(x, t) are written as 
E = —V(j) — c~^dA/dt and B = V x A, where the elec¬ 
trostatic potential and the vector potential A are de¬ 
termined by Poisson’s equation, 


VV(x,t) 



/a(x,V,t)d^V 


—47rtT, (2) 


and Ampere’s law, 

A nr 

v2A(x,t) =-jT, (3) 

c 

respectively. Here, the Coulomb (or transverse) gauge 
condition V ■ A = 0 is used and the current den¬ 
sity j = J2a^o.naUa = J2a I (or any 


vector field) is written as j = Jl + jr, where = 
-(47r)-iV/d3x'(V' • j)/|x - x'l and jt = (47r)-iV x 
(Vx/dVj/jx — x'l) represent the longitudinal (or ir- 
rotational) and transverse (or solenoidal) parts, respec¬ 
tively!^ Equations o, (ED, and © are the governing 
equations for the Boltzmann-Poisson-Ampere system. 

Suppose that fa, </, and A which satisfy Eqs. m- 
([3]) are given. Then, for the electromagnetic fields E = 
—V(/ — c~^dA/dt and B = V x A given from cj) and 
A, we consider the distribution function f^ which is the 
solution of the Vlasov equation. 


^ -f V- V-f ^ |E(x,t) + X B(x,t)| 
= 0 . 


A’ 


fa (X, V, t) 


(4) 


We also assume f^ to coincide instantaneously with fa at 
a given time to so that fY (x, v,to) = /a(x, v,to)- There¬ 
fore, equations obtained from Eqs. @ and © with fa 
replaced by fY also hold at to. In other words, fY, 
(j), and A satisfy the Vlasov-Poisson-Ampere system of 
equations at to. Note that the Vlasov-Poisson-Ampere 
system of equations can be derived from the variational 
principle using the action X defined by Eq. (1) in Ref. 29 
where its variation 61 associated with infinitesimal trans¬ 
formations of independent and dependent variables [see 
Eqs. (10) and (15) in Ref. 29] are explicitly shown in 
order to apply Noether’s theorem for obtaining conser¬ 
vation laws of energy and momentum. Now, let us use 
fY, 4>, and A to define the action integral I over a small 
time interval, to — h/2 < t < to + h/2, during which the 
Vlasov-Poisson-Ampere system of equations are approx¬ 
imately satisfied by fY, 4>, and A within the errors of 
order h. Then, neglecting the errors of higher order in h, 
the variation 61 can be written in the same form as in 
Eq. (15) of Ref. 29, 


61 =- 


ptQ-\-h/2 


dt / d^x 


/ to — hj2 


-dGo^(x,t) + V-dG^(x,t) 

(5) 


where 6Go and 6G^ are written as 

6G^{^,t) = SY 6tE-P^ ■6 xe, 

<5G'^(x, t) = Q/' 6tE - ■ 6 ke + S^6(t)-'SA- 6A. (6) 

Here, tE, 6 xe, 6(j), and 6A represent variations of t, x, 
(p, and A, respectively, while SY, S^, and 

'Ea are defined by 


='^J (x, v,t) -p ea</^ 


+ ^(-|V/>P + |Bp) 


= 


=J2j /„^(x,v,t) Qma|v|2-Pea</^ 

1 XdA dA \ 
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a 

+ ^ (|V<^I' - s') I + ^ [-(V</.)(V(A) 

+ ((VA) - (VA)^) • (VA)^ - -(VA)^l , 

C 

&t = -^V4., and ^ (b X I + il) , (7) 

respectively, where the superscript T represents the 
transpose of the tensor and I denotes the unit tensor. 
The field variable A which appears in Eq. © is intro¬ 
duced in Ref. 29 as the Lagrange undetermined multi¬ 
plier to derive the Coulomb gauge condition V • A = 0 
and it is shown from Eq. (8) in Ref. 29 that 


where tg is rewritten as t because tg is an arbitrarily 
chosen time. Equation (IT^ shows that 5Kcg represents 
effects of ICa in Eq. m on the conservation law. If ICa is 
given by the Coulomb collision term only, 6Kgo defined 
by Eq. (USD vanishes because the collision term conserves 
particles’ number, momentum, and energy. 

Energy and momentum balance equations can be de¬ 
rived from Eq. m using symmetries of the system under 
infinitesimal time and space translations as shown later. 
Before deriving them, we first consider the equation for 
the particle number density Ua = f fa(P'v which is ob¬ 
tained by taking the velocity-space integral of Eq. m as 

^ + V • {UaUa) = j ICad^yr. (14) 


V2A = 47rV-j. 


( 8 ) 


Suppose that the variations te, Sxe, dtp, and 5A. are 
such that SX = 0 holds for an arbitrary x-integral domain 
in Eq. Then, taking the small time interval limit 
h —>• -1-0 in Eq. m, we find that the conservation law, 


d 




+ V ■6G^{x,to) = 0, 


(9) 


J t=to 


should be satisfied. This is the so-called “Noether’s the¬ 
orem. Recalling that fY (x, v, to) = /a(x, v,to) and com¬ 
paring Eqs. o and ® at t = to: we have 


'dfa (x,V,t)' 


'dfa{x,V,t)' 

9t 

t=to 

dt 


/Ca(x, v,to)- 


( 10 ) 

We now define £c, Pc, Qc, and He from P^, Q)f, 
and n^, respectively, by replacing with fa in Eq. ©■ 
Correspondingly, 5Gg and (5G are defined from SGg and 
by replacing £1^, P)f, Q)f, and 11^ with £c, Pc, 
Qc, and He, respectively, in Eq. (0. These definitions 
immediately yield SG^{x,tg) = (5G(x, to) and 


dSGg (x, t) 


'dSGg{x,t))' 

dt 

t = to 

-1 

_1 


6KGo{x,to), 

( 11 ) 


where Eq. m is used and SKgo is defined by 


(5Ago = Ksc dts — Kpe • 6xe, 

Ksc = X) / , 

Kpe = ^ ICa {maV + ^a) . (12) 

a 

Substituting Eq. (HU into Eq. we find that the 
conservation law is modified for the Boltzmann-Poisson- 
Ampere system as 




SGg{x, t) -b V ■ ^G(x, t) = SKgo, 


(13) 


We hereafter assume that Ca f ICad^v = 0, which 
means that the source terms ICa conserve electric charge 
even if J ICad^v 0 for each species a. This seems a rea¬ 
sonable assumption in consistency with Eqs. m and m 
in which no external source terms are included. Then, 
multiplying Eq. (I14|) with the electric charge Cq and per¬ 
forming the summation over species result in the charge 
conservation law. 


| + V.J = 0. (15) 

We also find from J2a f l^ad^v = 0 that, in Eq. (fT^ . 
the terms including (j) and A vanish and make no contri¬ 
bution to Afc and Kp^. As seen from Eqs. ®, and 
(HU), we can put A = d(j)/dt which is also used in Ref. 29 
to derive energy and momentum conservation laws for 
the Vlasov-Poisson-Ampere system. 

We now note that the action integral is invariant, 
namely, SI = 0 under the infinitesimal translations in 
space and time represented by dtp = cq, Sxa = Sxe = e, 
Sva = 0, Sej) = 0, and SA = 0, where eg and e are constant 
in time and space. These invariance properties hold be¬ 
cause the integrands in the action integral X depend on 
(x,t) only through variational variables [see Eq. (1) in 
Ref. 29]. Using the time translational symmetry, Eq. (fT51) 
reduces to the canonical energy balance equation, 

d£ 

^ + V • Q, = Kea, (16) 

where the canonical energy density and flux {£c, Qc) 
are given by replacing f^ with fa in the definitions of 
{£Y, Q^) in Eq. ([71). In the same way as in Ref. 29, we 
use the kinetic energy density and flux, {£p, Qp), defined 

by 

/a(x,V,t)ima|vp, 

Qp=X//a(x,V,t)ima|vpV, (17) 
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to modify Eq. dm) into more familiar forms. Then, the 
energy balance equation is finally written as 

a , |Ei|2 + |Bp 

dt V ^ + 8^ 

+ v.(q, + ^e><b-T|e,) 

d , |Ei|2 + 2Ei.ET + |B|2 

“ dt Stt 

+ V • (Qp + ^E X B + 

= isTgc = /Ca^malvp, (18) 

a 

where E^ = —V(j> and Et’ = —c~^dA/dt are the longi¬ 
tudinal and transverse parts of the electric field, respec¬ 
tively. 

Next, from the space translational symmetry and 
Eq. ()13p . we obtain the canonical momentum balance 
equation, 

dP 

^ + V • n, = Kpc, (19) 

where the canonical momentum density and tensor 
{Pc^c ) are given by replacing with fa in the def¬ 
initions of (P)f , Tla ) in Eq. ©■ Eurthermore, in the 
same way as in Ref. 29, the invariance of I under the in¬ 
finitesimal rotation is shown to give the equation for the 
angular momentum, which is used to modify Eq. m 
into the momentum balance equation, 


^(Pp + P/)+V-(np+n/) = Kpa = Y, j ICamaV. 

( 20 ) 

Here, the particle parts (Pp,!!^) of the momentum den¬ 
sity and the pressure tensor are defined by 

Pp = ^ / /a(x, V,t)maV, 

a 

rip = d^v /a(x,V,t)maVV, (21) 


and the field parts (P/,!!^:) are given by 


P/ = 

1/ 


Er X B 


47rc 


1 

87r' 


Uf = ^(lErP 


■ 2 EL-ET + |Bni 


— — (E^Ei 4-EiEx 4 -EtE^- b BB). (22) 

47r 


Equations (fT^ and (j20() take physically familiar forms 
of energy and momentum balance equations including ex¬ 
ternal source terms. As mentioned earlier, if fCa is given 
by the Coulomb collision term only, Kgc and Kpc van¬ 
ish so that the energy and momentum balance equations 
for the Boltzmann-Poisson-Ampere system take the same 
forms as those for the Vlasov-Poisson-Ampere system<^ 


III. GYROKINETIC BOLTZMANN-POISSON-AMPERE 
SYSTEM 


Let us start from the gyrokinetic Boltzmann equation 
written as 


m 


dZa 

dt 



Fa{Z,t) 


= E^aJ^a,E,](Z,t)4-5„(Z,t), 

b 


(23) 


where Fa{Z,t) is the gyrocenter distribution function for 
species a, Cai,[Fa, Fb]{Z,t) represents the rate of change 
in Fa{Z,t) due to Coulomb collisions between particle 
species a and b, and iSo(Z,t) denotes other parts includ¬ 
ing external particle, momentum, and/or energy sources 
if any. The gyrocenter coordinates are written as Za = 
(Xa,t7a,Aia,Ca), where Xa, Ua, fJ-a, and fa represent the 
gyrocenter position, parallel velocity, magnetic moment, 
and gyrophase angle, respectively. Appendix A shows 
the relation of the gyrocenter coordinates to the particle 
coordinates in detail. The perturbation expansion pa¬ 
rameter in the gyrokinetic theory is denoted by 6 which 
represents the ratio of the gyroradius p to the macro¬ 
scopic scale length L of the background field. It is shown 
in Appendix B how the collision operator (7®^ [Fa , Fb] for 
the gyrocenter distribution functions Fa and Ft is given 
from the collision operator Cabifa, fb] for the particle dis¬ 
tribution functions fa and fb- 

The deviation of each distribution function from the 
local Maxwellian is regarded as of 0{S), and accord¬ 
ingly the collision term is considered to be of 0{5). 
We assume that the source term Sa is of 0(6'^) so that 
its effect appears only in the transport time scale. We 
also assume that J^a^aJ dU J dp J df DaSa{Z,t) = 0 
in order to prevent the source term from affecting the 
charge conservation laws [see Eq. (15(1)) ]. Here, Da de¬ 
notes the Jacobian for the gyrocenter coordinates. Da = 
det[a(xa, Va)/9(Xa, Ua,fa,Pa)], where (Xa, Va) represent 
the particle coordinates consisting of the particle’s posi¬ 
tion and velocity vectors. 

We treat toroidal systems, for which the equilibrium 
magnetic field is given in the axisymmetric form as 


Bo = V X Ao = /VC + VC X Vx, (24) 


where I and x are constant on toroidal flux surfaces la¬ 
beled by an arbitrary radial coordinate s and C is the 
toroidal angle. We note that I and x represent the co¬ 
variant toroidal component of the equilibrium field Bg 
and the poloidal magnetic flux divided by 27r, respec¬ 
tively. The equilibrium field Bq is allowed to be depen¬ 
dent on time. Then, following Ref. 18, the gyrocenter 
motion equations are written as 

§l = {Z.,aJ + (Z.,X4.£^^, (25) 

where the gyrocenter Hamiltonian, which is independent 
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of ^a, is defined by 


1 9 

Ha = -^rUaUa + ^J.aBo + Ca'ifa, 


(26) 


and A* is given by A* = Ao(Xa,f) + 
{maC/ea)Uah('X.a,t). Using the nonvanishing com¬ 
ponents of the Poisson brackets for pairs of the 
gyrocenter coordinates given by 


{X„XJ = —^bxl, {Xa,Ua}= 


eaB* 


rUaB* 




(27) 


the gyrocenter motion equations in Eq. (I25p are rewritten 
as 

dX, 1 


dt 




m, 


a dUa 


^all L 


+ Cb X ( 

Cn C Ot 


dUg 

dt 


b: 


m„B*u 


( 1 dK* 

Bq + Ba ( V'I'a + ~ 


d^.a 

dt 


= 0 , 


and 


d£.a ^ el d'ita 

dt nriaC d^a 


, (28) 


(29) 


(30) 


(31) 


Here, ila = eaBol{niac), b = Bq/Bq, H*|| = B* • b, and 
B* = V X A*. The field variable d'a is defined by 

Vk, = (V;„(Z„,t))^^ + (|Ai(X, 

where the field variable 'i/'a is defined in terms of the 
electrostatic potential (j) and the perturbation part of the 
vector potential Ai as 

ijjai'^a.t) = (l){yia+Pa,t) “ -Vao(Za,t) • Ai (Xa -f : *) ■ 

(33) 

The gyroradius vector is given by = b(XQ,t) x 
Vao(Za, t)/na(Xa, t) and the zeroth-order parti¬ 
cle velocity v^o is written in terms of the gyro¬ 
center coordinates as Voo(Zo,t) = C/ob(Xa,t) — 
[2poBo(Xa)/ma]^/^[sm,faei(Xa,t) -f C0S.foe2(Xa,t)], 
where the unit vectors (ei,e 2 ,b) form a right-handed 
orthogonal system. The gyrophase-average and 
gyrophase-dependent parts of an arbitrary periodic 
function Q{^a) of the gyrophase ^a are written as 

^Q(ea) and Q = Q-{Q)^^, (34) 


respectively. In the gyrocenter motion equations, ef¬ 
fects of the time-dependent background magnetic field 
and those of the fluctuating electromagnetic fields ap¬ 
pear through dA*/dt and 4>a, respectively. It should 
be noted that dZa/dt on the left-hand side of Eq. (l23l) 
is regarded as a function of (Z,t) which is given by the 
right-hand side of Eq. (|25|l . 

We find from Eqs. (IMl)-(IX4l) in Appendix A and 
Eq. (IBII) in Appendix B that the gyrophase-dependent 
part of the right-hand side of Eq. (1^ appears from 
and it is of 0{S). Using Hq = 0{S~^), the gyrophase- 
dependent part of the left-hand side of Eq. (1^ is written 
as rtadFa/d^ to the lowest order in S. Then, it is con¬ 
cluded that Fa = 0{S^). Taking the gyrophase average 
of Eq. (1231) , we obtain 


m 


dZa 

dt 


A 

ffZ 


FaiZ,t) 


= J2{C^ab{Fa,F,]iZ,t))^+SaiZ,t), (35) 

a 

where Fa{Z,t) and SaiZ,t) are both regarded as inde¬ 
pendent of the gyrophase ^ and (■••){ are omitted from 
them for simplicity. It is seen from Eq. m that ef¬ 
fects of Fa = 0(i5^) on {C^i,[Fa, Fb])^ in the right-hand 
side of Eq. (1551) are estimated as of 0{6^). Here and 
hereafter, we neglect Fa = 0{6^) in both sides of the gy- 
rokinetic Boltzmann equation given by Eq. ()35|) . Even 
so, its moment equations can correctly include the colli- 
sional transport fluxes of particles, energy, and toroidal 
momentum up to the leading order, that is 0(5^), as con¬ 
firmed later. In Appendix C, Eq. m combined with 
Eqs. (IC2I) . (ICIOI) . (IC15|) . and (ICI6I) presents the approx¬ 
imate gyrokinetic collision operator, which has favorable 
conservation properties and correctly describes collisional 
transport of energy and toroidal angular momentum. 

The gyrokinetic Poisson equation and the gyrokinetic 
Ampere’s law are written asi^ 

V2(/,(x,t) = f d^Z Da{Z,t)S^{X + Pa - ^) 


FaiZ,t) 


Bo dp, 


and 


/L'jr 

V (Ao-|-Ai) =-(jG)T5 

c 


(36) 


(37) 


respectively, where (jG)T is the transverse part of the 
gyrokinetic current density jc defined by 

= / d^ZDa(Z)S^lX + Pa(Z)-x] 


X [^Fa(Z,t) 

Ba'4^a dFa 
Bo dp 


v,o(Z)-^Ai(X + pJZ),t) 

rriaC 


Vao(Z) . 




(38) 
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It should also be noted that the Coulomb gauge condi¬ 
tions V • Aq = 0 and V • Ai = 0 for the equilibrium 
and perturbation parts of the vector potential are used 
here. The equilibrium vector potential Aq is given by 
Ao = —X^C+VC X V77, where 77 = ri{R, Z) is the solution 
of A*77 = • {R~‘^S/r]) = I. In Ref. 18, additional gov¬ 

erning equations are derived in order to self-consistently 
determine / and x for the time-dependent axisymmetric 
background field Bq = V x Aq = /VC + VC x Vx. They 
are given by 


I = 


dO 

27r 


47r 


Me-h - Ric 


(39) 


and 


A,x = 


47r 


[(j(gc))j, + V X M] 


+ 

dx 


- V X B 


■/?2VC 

(40) 


where the toroidal-angle average is represented by • • • = 
(27r)~^ j> • • • dC, the poloidal angle is denoted by 0, and 

Ae = / — (47r/c)Me — (i?fo'=))e + Ri(- Here, the covariant 
toroidal component of an arbitrary vector V is written 
as Ve- The turbulent part of the magnetic field is given 
by Bi = V X Ai and is defined by V x B*^sc) = 

(47r/c)(j^s‘^^)T. The magnetization M can be obtained 
from the turbulent fields and the distribution functions 
for all species using Eqs. (41)-(43) in Ref. 18. Thus, 
Eqs. (1551) . (1551) . (1571) . (15^ . and (150)) constitute the closed 
system of governing equations which determine Fa, (j>, 
Ai, I, and x- 

For the gyrocenter coordinates which have Poisson 
brackets given by Eq. (EZl), the Jacobian is given by 
Da = B*^^/ma- It is important to note that the Jacobian 
Da satisfies the gyrocenter phase-space conservation law. 


dDa 

dt 


A 



(41) 


Then, using Eq. m, the gyrokinetic Boltzmann equa¬ 
tion in Eq. dSSl) can be rewritten as 

I {DM + A . = DM. (42) 


where ICa is the gyrophase-independent function given by 
the right-hand side of Eq. (1551) , 


ICaiZ,t) =Y,{C^dFa,F,]iZ,t))^ + SaiZ,t). (43) 

a 


We hereafter derive conservation laws for the gyroki¬ 
netic Boltzmann-Poisson-Ampere system of equations 
following the procedures similar to those shown in Sec. II. 
For that purpose, suppose that Fa, 4>, Ai, I, and x satisfy 
Eqs. (1551) . (1^ . (1571) . (1551) . and (1501) . Then, we consider 


the gyrocenter distribution function F^ which obeys the 
gyrokinetic Vlasov equation. 




dZa 

dt 



fY = 0, 


(44) 


where dZa/dt is evaluated by using the above-mentioned 
fields (()),Ai,/,x) obtained from the solution of the 
gyrokinetic Boltzmann-Poisson-Ampere system of equa¬ 
tions. Here, it should be noted that, if the distribu¬ 
tion functions Fa and FY , which are given as the so¬ 
lutions of Eqs. (I35|) and (1551) . respectively, are initially 
gyrophase-independent, they are gyrophase-independent 
at any time. Besides, FY is assumed to coincide instanta¬ 
neously with Fa at a given time to. Therefore, Eqs. dMl), 
§3), OMl), and (1551) are all satisfied at that moment even 
if Fa is replaced with FY in these equations. Thus, the 
gyrokinetic Vlasov-Poisson-Ampere system of equations 
are instantaneously satisfied by (E))^, 0, Ai,/, x) at t = 
to. In Ref. 18, the action integral I is defined to derive 
all the governing equations for the gyrokinetic Vlasov- 
Poisson-Ampere system based on the variational princi¬ 
ple, and its variation 61 associated with the infinitesimal 
variable transformations are given to obtain conservation 
laws from Noether’s theorem. Here, the action integral X 
can be expressed in terms of {FY, (j), Ai, /, x) over a small 
time interval, to — h/2 < t < to + h/2, during which the 
gyrokinetic Vlasov-Poisson-Ampere system of equations 
are approximately satisfied by them within the errors of 
order h. Then, neglecting the errors of higher order in 
h, we can write the variation SI in the same form as in 
Eq. (77) of Ref. 18, 


^JG^(X,t) + V-JG^(X,t) 

(45) ■ 

with the functions SGq and defined by 

dG^iX,t)=£Y StE-F^ -S^E, 

JG^(X, t) = QY StE - nY ■ d^E + S^d^- -Eai ■ SAi 
~i:^o-SAo + S^6x + ST^, (46) 

where £Y and are defined by 


SI=- 


lO-l-D/Z 


dt / d^X 


tr)—hl2 


E 


J dU J d^iJ d^ DaFYHa 


+ ^(-|V0P + |Bo + Bi|2), 

Y,JduJd^lJd^ DaFY {maUh + YlAo) , (47) 


and definitions of other variables nr, s^, Sai, sXo 
S^, and (5T'^ are shown in Eq. (79) of Ref. 18. The super¬ 
script V in the variables (£Y,PYy' implies that they 
are defined using the distribution function FY instead of 
Fa- 

As explained in Ref. 18, the integral domain of Eq. (1551) 
is not an arbitrary local one in the X-space but it can be 
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local only in the radial direction in order for Eq. (1451) to 
be valid. Then, if the variations 6tE,SxE, ■'' in Eq. dUl) 
are such that 61 = 0 holds for a spatiotemporal in¬ 
tegral domain defined by [to — h/2,to + h/2] x [si,S 2 ] 
where [si,S 2 ] represents an arbitrary spatial volume re¬ 
gion sandwiched between two flux surfaces labeled by si 
and S 2 , then the conservation law is derived as 


^Go^(X,t)-KV-<5G^(X,to) 


t=to 


pGUX.t)\ + j;^JV'{SG''-Vs)) 


= 0 . 


t=to 

(48) 


This is Noether’s theorem for the gyrokinetic Vlasov- 
Poisson-Ampere system. In Eq. (H51) . V' = dVjds rep¬ 
resents the derivative of V(s, t) with respect to s and 
V{s,t) denotes the volume enclosed by the flux surface 
with the label s at the time t. 

Using {Z,to) = Fa{Z,to) and comparing Eq. (l35ll 
with Eq. (IT41) . we find 


\dFnZ,t)^ 


\dFa{Z,t)-\ 

dt 

t = to 

dt 


ICa{Z,to), (49) 


their shapes and the grid of the flux coordinates moves. 
Then, Eq. (1551) is rewritten as 

- (U' (^Go))+^ (U' {{SG - SGou,) ■ Vs)) = U' {SKgo) , 

(54) 

where Ug ■ Vs represents the radial velocity of the 
flux surface labeled by s and Ug is defined by Ug = 
dx{s,9X,t)/dt with the flux coordinates (s, 0, C) [see 
Eq. (2.35) in Ref. 31]. In Sec. V, gyrokinetic energy 
and toroidal angular momentum balance equations are 
derived from Eq. (l54)) . 


IV. EQUATIONS FOR GYROCENTER DENSITIES AND 
POLARIZATION 


In this section, we take the velocity-space integral of 
the gyrokinetic Boltzmann equation in Eq. (I35p to con¬ 
sider the particle transport before treating the energy 
and toroidal angular momentum transport in Sec. V. We 
define the gyrocenter density by 

n(s=)(X,t) = JdU JdtiJd^ DaF,, (55) 


where /C is defined by Eq. (PI). Let us also define 6Go 
and SG from SGq and <50^^ by replacing F^ with Fa. 
Then, we have G^(X, to) = G(X,to) and 


'ddG^{X,ty 


r5(5Go(X,t))1 

dt 

t=ta 

dt 


where 


SKgoO^, to), 
(50) 


SKgo = Ksc6tE — Kpc • 

Ksc = '^ea f dU f dfi f d^ DaK-aHa 


Kp, = ^ e< 


f dU j d^i f d^ DalCaPa- (51) 


a 

Here, p° denotes the canonical momentum for species a 
defined by 


= ^A: = ^Ao + maUh. (52) 

c c 

Substituting Eq. (ISOl) into Eq. (H51) and rewriting the ar¬ 
bitrarily chosen time to as t, we obtain the conservation 
law for the gyrokinetic Boltzmann-Poisson-Ampere sys¬ 
tem, 

(^^^SGo{X,t) + V-SG{X,t)'^ 

= {SKgo) , (53) 

where (SKgo) represents effects of the collision and 
source terms on the conservation law. Under the nonsta¬ 
tionary background field Bq, flux surfaces may change 


and the gyrocenter flux by 

= JdU JdfiJ diDaFaV^r\ (56) 


where represents the gyrocenter fluid velocity and 
the gyrocenter drift velocity = dXa/dt is given by 
evaluating the right-hand side of Eq. (I28p at (X, [/,/i). 
Then, integrating the gyrokinetic Boltzmann equation, 
Eq. (1421) . with respect to the gyrocenter velocity-space 
coordinates (C/,/r,^) and using Eq. (HSl) . we obtain 

+ V • = J dU J dfl J d^ DaSa. 

(57) 

Using the approximate collision operator given by 
Eq. (ICll) in Appendix C, the particle flux F^ due to 
collisions and finite gyroradii is defined by Eq. (IC6I) with 
putting AP(z) = 1. [If the collision operator given by 
Eq. (IBII) in Appendix B is employed, F^ is defined by 
Eq. (lB8l) .1 

As shown in Ref. 18, the gyrokinetic Poisson equation 
in Eq. (1551) is rewritten as 

= (58) 

n. ^ 


where Ep = — V^, V = d/dX, and p(P°6 represents the 
polarization density defined by 




(- 1 )- 


n—1 

^n-1/ 


n\ 


d^-^{DaF*p^Pa,, ■ • 


j dU j dp j d^ 

(59) 





















Here, pai denotes the ith Cartesian component of 
Pa = b(X,t) X Vao(Z,t)/Ha(X,t), and = Fa + 

{ea'il^a/Bo){dFa/dp). 

As mentioned before Eq. (Oil) in Sec. Ill, 
J2a^af dU J dpj DaSa{'Zi,t) = 0 is assumed. Then, 
using Eq. (113, we can obtain the charge conservation 
law. 


- (y 

dt 




+ V-(j(s=)+f)=0, (60) 


where the current density due to the gyrocenter drift and 
that due to the collisional particle transport are given by 
jtec) = = Ea eany uy and = E, e.r^, 

respectively. Note that the magnetization current is 
solenoidal and, accordingly, it does not contribute to the 
charge conservation law in Eq. O- Equation (1551) is 
substituted into Eq. (l60l) to show 


dgc) I :C _ 

Jl Jl 


47T 


py ^ ) , (61) 


where the subscript L is used to represent the longitu¬ 
dinal part of the vector variable. Then, using Eqs. 
dini) and dUl), we find that the useful formula, 


d 




(gc) 




= V- 


E 


Con 


(gc) 


+ (jy •V-4)-(A(V-jg)) 


dt I 47r 


>(pol) 


\dt 


47r ^ 


(pol) 




VA 


(62) 


holds for any function A(X, t). The relation in Eq. (1551) 
is used in Sec. V.B to derive Eq. dZl). 


V. GYROKINETIC ENERGY AND TOROIDAL 
ANGULAR MOMENTUM BALANCE EQUATIONS 


and they satisfy Eq. (1541) which, in the same manner as 
in Ref. 18, leads to the energy balance equation, 

^ (E' (f)) + ^ {V ((Q: + - fu,) • Vs)) = E' {Ksa) 

(63) 

Here, the energy density £ is defined by 


£ = £c + V ■ [ 


= J2jdUjdpjd^ DaF, 
+ T(|V0|" + |B„ + B,|") 


mg 

2 


VaO - 


-Ai 


(64) 


and the energy fluxes Q* and QJj are given by 




^ dU dp d^ DgFa 




OAq 

dt 


X ( -ph + ^^^(vte'=))_L - Na 


Bo 

1 9(Ao + Ai) 
47r dt 
i9Ao 


1 


x(b, + b,) + 3T(x^ 


Q- 


dt 




and 


Qfl = Qi 


dt 


(65) 


( 66 ) 


respectively, where N^, Q/j, and are defined by 
Eqs. (43), (85), and (88) in Ref. 18, respectively. 

Recalling Eq. m and using Eq. (IC5I) into which we 
substitute A®(Z) = F[a{Z) and the approximate collision 
operator given by Eq. (ED in Appendix C, we rewrite 
the right-hand side of Eq. (1551) as 


In this section, energy and toroidal angular momentum 
balance equations for gyrokinetic systems including col¬ 
lisional processes are derived by using the results shown 
in Secs. HI-IV. 


A. Energy balance equation 

The variation 51 of the action given in Sec. HI van¬ 
ishes under the infinitesimal time translation represented 
by StE = e where e is an infinitesimally small constant. 
Here, all other infinitesimal variations 6xe, 54), ■■■ are 
regarded as zero. Then, 5Go and (5G are determined 
by these conditions for the infinitesimal time translation 


P'(A£e) = --(E'(Q^-Vs)) 


V'Y,( dU dp d^DgSaHa 


(67) 


where the energy flux due to collisions and finite gyro- 
radii is defined by taking the summation of Eq. (IC6|) over 
species a with putting Ag(z) = ^rngV^^^ -I- poaBo{xa) + 
ea4>{^a)- The derivation of Eq. (1571) requires Eq. (IC7I) 
which is satisfied by the collision operator in the form of 
Eg. dCip with appropriately choosing Axi^\ and 

Apoa as described in Appendix C. 
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Substituting Eq. (EH) into Eq. (|63)) , the energy balance 
equation is rewritten as 


(^)) + ^(E'((Q-£:u,)-Vs)) 
= V'Y,( [ dU [ dfl [ DaSaH, 


( 68 ) 


where the energy flux Q is given by 

Q = Q: + Qfl + Q'^- (69) 

The right-hand side of Eq. (|68)) represents the external 
energy source. It is confirmed later in Sec. VLB that 
the ensemble average of (Q • Vs) coincides with the well- 
known expression of the radial energy transport to the 
lowest order in the ^-expansion. 


B. Toroidal angular momentum balance equation 

The toroidal angular momentum balance equation is 
derived from the fact that 5X = Q under the infinitesimal 
toroidal rotation represented by dxs = ee^(X). Here, e 
is again an infinitesimally small constant, and e(^(X) is 
defined by e^(X) = d'X{R,z,()/dC = where the 

right-handed cylindrical spatial coordinates {R, z,() are 
used. We also define z by z = i?VC x Vi? which repre¬ 
sents the unit vector in the ^-direction. Then, if putting 
the origin of the position vector X at (i?, z) = (0,0), we 
have e^(X) = X x z. Under the infinitesimal toroidal 
rotation, the variations of the vector variables are given 
as (5Ai = eAi x z and (5Ao = cAq x z although the other 
variations dtE, 54>, •••, are all regarded as zero. Then, 
using these variations of the variables associated with the 
infinitesimal toroidal rotation, the canonical momentum 
balance equation is derived from Eq. (1531) as 


the variation of the canonical toroidal angular momen¬ 
tum due to collisions and external sources is given by 


K 


Pci 


= K 


Pc • 


ec = ^IdC/y dC DalCa{pl)c. 


(73) 


We follow the same procedures as shown in Sec. V.B 
of Ref. 18 and use Eqs. dZoJ-dzl and Eq. (I62|) with 
A = Ao(^ = —X to write the toroidal angular momen¬ 
tum balance equation as 


dt 


P|ic--(P 


>(pol) 


47r 


• VA, 


oc 


j__a 

^V'Ws 


1 


W { Hf^ +^RC-^ Mic(V X Bi) • Vs) 




dA, 


oc 


dt 


j(pol) 




(Kpcc) + ■ (^ocJl)) 


(74) 


where 


^iic = 

a 


dU 


dfi 


d^ DaFaTUaUbi, 


a 

= Pflc - Vs. 


dC DaFamaUbcvi^^^ ■ Vs, 

(75) 


a(P,-ec)\ , 1 d 


dt 


V' ds 


[U' (Vs • (n, • ec 


+ (^Ai X Ai -|- X Aq) • z -L P i{(^))] 
= {Kpcc) ■ 


(70) 


Here, the density of the canonical toroidal angular mo¬ 
mentum is defined by 


Pc • = 


] J dU J dp J dC DMp:)c, (71) 


with the toroidal component of the canonical momentum 
denoted by 


{p:)c = ^A:^ = ^Aoi + maUbc, (72) 

where b^ = I/Bq represents the covariant toroidal com¬ 
ponent of b = Bo/Ho- Definitions of S^i, S^iq, and P^^ 
on the left-hand side of Eq. dZOl) are given by Eqs. (79) 
and (100) in Ref. 18. On the right-hand side of Eq. (1701) . 


Using Eqs. (173]) and (IC5p in Appendix C with putting 
Aa = {Pa)i^ the right-hand side of Eq. (174)) is rewritten 
as 


(A>o<) + 1 <V . (Ao^jJ)) 

‘ " U((jc +ijs 


V ds L 


■Vs 


dU dp d^ DaSarUaUbi ) , (76) 


where is defined by taking the summation of 

Eq. (IC6p over species a with putting A0(z) = 

(ea/c)[Aof(x) -I- Ai,;(x)] -I- rriaVi. In deriving Eq. ( [76l), 
dU J dp J d^ DaSaiXjA) = 0 and Eq. (10121) are 
used. The approximate collision operator which satisfies 
Eq. ()C12D is presented in Appendix C. 

Substituting Eq. (1701) into Eq. (170)) . the toroidal angu- 
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lar momentum balance equation is rewritten as 




dt 


V ds 


V {+ u 


RC 


4tt 

C\s 


'VA, 


oc 


- s; 


X (A,((V X Bi) ■ Vs) - — {EitEi + 

47r 


+ 


1 


47rc \ d( 


dX 


A1 


1 ! dA 


oc 


dt 


>(pol) 


+ 




•Vs 


%- 


1 


,(pol) , 




oc 


47r ) 

(Us • Vs) 


= '^(^J du J d^i J DaSamallb^'^ , (77) 


We assume that the ensemble average (Q)ens of 
any variable Q considered here has a slow temporal 
variation subject to the so-called transport ordering, 
d\n{Q)ens/dt = 0{6 ^vt/L), and that it has a gradient 
scale length L which is on the same order as gradient 
scale lengths of the equilibrium field and pressure pro¬ 
files. We also impose the constraint of axisymmetry on 
{Q)ens that is written as d{Q)ens/dC = 0 even though Q 
itself is not axisymmetric. On the other hand, the turbu¬ 
lent part Q of Q is assumed to vary with a characteristic 
frequency oj = 0{vtIL) and have gradient scale lengths 
L and p in the directions parallel and perpendicular to 
the equilibrium magnetic field Bq, respectively. 

The ensemble-averaged part {Fa)ens of the distribution 
function Fa for species a consists of the local Maxwellian 
part and the deviation from it, 


where the right-hand side represents the external source 
of the toroidal angular momentum and 

= (jC + . Vs (78) 

is the radial flux of the toroidal angular momentum due 
to collisions and finite gyroradii. In Sec. VI.C, we derive 
the ensemble-averaged toroidal angular momentum bal¬ 
ance equation from Eq. (EZl) in order to confirm that it is 
consistent with the conventional result up to the second 
order in S. 


VI. ENSEMBLE-AVERAGED BALANCE EQUATIONS 
FOR PARTICLES, ENERGY, AND TOROIDAL ANGULAR 
MOMENTUM 

In this section, the particle, energy and toroidal an¬ 
gular momentum balance equations derived in Secs. IV 
and V are ensemble-averaged for the purpose of verify¬ 
ing their consistency with those obtained by conventional 
recursive formulationsj ^°i^^~ — In the same way as shown 
in Sec. VI of Ref. 18, we divide an arbitrary physical 
variable Q into the average and turbulent parts as 

Q = (Q)ens + Q, (79) 

where (• • • )ens represents the ensemble average, and we 
immediately find (Q)ens = 0. We identify the zeroth 
fields Ao and Bq with the ensemble-averaged parts to 
write Ao = (A)ens, Ai = A, Bq = (B)ens, and Bi = B. 
Regarding the electrostatic potential (j), it is written as 
the sum of the average and fluctuation parts, ^(x, t) = 
((()(x,t))ens+0(x,t). Here, assuming that {(j){x,t))ens 0, 

the background E x B flow is retained and its veloc¬ 
ity is regarded as 0{5vt), where S and vt represent the 
drift ordering parameter and the thermal velocity, respec¬ 
tively. Then, using Eq. ([33]), we have V'o = (V’a)ens +'0a, 
where 

0a = 0 - — • A. (80) 
c 


(Fa) ens — ^aM + (Fal) ens ■ (81) 


The local Maxwellian distribution function is written as 
FaM = nao[ma/{2TTTao)f/‘^exp[-{maU‘^/2 + pBo)/Tao] 
where the equilibrium density Uao and temperature Tao 
are regarded as uniform on flux surfaces. The first-order 
ensemble-averaged distribution function {Fai)ens is de¬ 
termined by the drift kinetic equation, which can be de¬ 
rived by substituting Eq. (ISTj) into the ensemble average 
of Eq. (1^ . The derived equation agrees, to 0{S), with 
the well-known linearized drift kinetic equation, on which 
the neoclassical transport theory is based^ii^^ 

The fluctuation part Fa is written as 


Ta — —-faM --h Ra- 

J- n. 


(82) 


Substituting Eq. (l82|) into the fluctuation part of the gy- 
rokinetic equation in Eq. (1231) yields 


dhg 

dt 

= FaM 


+ {ha, Ha} 

Ga di^pa)^ x(gc) Sa . , 


TaO 

irUaU'^ + p.Bo 5 


F 


aO 


- - ) V In Tao 




(83) 


where Cg represents the linear collision term defined by 


Cf'(X) = Y,{[Cl,{ha{x - Pa),F,M) 

b 

+ F!ai,iFaM, hb[x - Pb))]^=X+pJi^- (84) 

Equation (1831) is valid to the lowest order in 5 and agrees 
with the conventional gyrokinetic equation for the nona- 
diabatic part ha of the perturbed distribution function 
derived from using the WKB representation! On the 
right-hand side of Eq. (155)) . the turbulent part of 
the gyrocenter drift velocity = {Xq, Ha} 

is written as 


X V(0,(X + Pa, t))^ + 0{S^). (85) 


{'4^a)ens — (0)ens; 
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It is shown by using Eq. (1821) and the WKB representa¬ 
tion that, to the lowest order in d, the turbulent parts 
of Eqs. (|36ll and (1^ agree with the gyrokinetic Poisson 
equation and the gyrokinetic Ampere’s law derived by 
conventional recursive formulations^i^S 


A. Ensemble-averaged particle balance equation 


Taking the ensemble average of Eq. (ElD and subse¬ 
quently its flux surface average, we obtain 



where {n^a^^)ens = nao -I-0((5), and ((•••)) represents 
a double average over the flux surface and the ensem¬ 
ble. Here, nao is the equilibrium density which is a flux- 
surface function and characterizes the Maxwellian distri¬ 
bution function FaM- On the right-hand side, the source 
term Sa is regarded as of 0{S‘^) as well as all other terms 
in Eq. (1861) . and it is assumed to have no turbulent com¬ 
ponent so that Sa = {Sa)ens- 

It is shown in Ref. 18 that the radial gyrocenter particle 
flux is given by 

(r(gc)). = . Vs)) = + (rA)^ (s?) 

where the nonturbulent part and the turbulence- 

driven part are written as 

(T^^y = QdU JdfiJ de Da{Fa)ens{yri«^^)ens ' 

+ nao {mens Xh)-Vs^+ 0 {S^), (88) 

and 


{T^y = QduJd^iJ d^ DaiFavi^mns ' 



/"G 

[ d^ DahaiVlpa X b) 

•Vs' 

)> 

+ 0(^3), 




(89) 

respectively. 

On 

the right-hand 

side 

of 

Eq. (Ell), PCGL 

represents the hrst-order 

part 


of the pressure tensor in the Chew-Goldberger- 
Low (CGL) formal defined by = 

JdUjdfiJd^ Da{Fal)ens [rUaU^hh + flBo {1 - hh)] , 
and the ensemble-averaged electric field is given by 
(E)ens = — V((/))ens “ c~^dAo/dt. Thus, Eq. (1551) ex¬ 
presses the neoclassical radial particle flux and the radial 
E X B drift which are well-known by the neoclassical 


transport theory^i We also hnd that Eq. (I89p agrees 
with the turbulent radial particle flux derived from the 
conventional gyrokinetic theory based on the WKB 
formalismi^ 

The radial classical particle flux is given by 


{T^y = {(T^-Vs)) 

= E JduldfiJd^ Da[{v X b) • Vs] 

X [C^abihmbM) + C^yFaMjbl)]) + 0{6^) 

= /m^{Faixh)-Vs)+0{d^), (90) 

\ eaBo / 


where Fai = f d^'vma'vC^ is the collisional friction force. 
It is well-known that the classical transport equation re¬ 
lating (r^)® to the gradient forces is immediately de¬ 
rived from Eq. (IMl) because the first-order gyrophase- 
dependent part of the particle distribution function in 
Eq. (|M1) is expressed in terms of the gradient of the 
background Maxwellian distribution function as fai = 
~Pa ■ ^FaM with the gradient operator V taken for the 
fixed energy variable e = ^niaV^ + e{(j>)ens- 

In the same manner as in deriving Eq. (I54|) from 
Eq. (1551) . the ensemble-averaged particle transport equa¬ 
tion can be obtained from Eq. (l86l) as 


^ {V'riao) + ^ (P' [(rj® - n,o(w • Vs)]) 


dU dp, d^ DaSa 


(91) 


where the total radial particle flux is given by the sum of 
the classical, neoclassical, and turbulent parts as 

(P,)® = (Pte^))® + (pC)® 

= (P^A)® + (pA)® + (pC)«. (92) 

As shown above, the well-known expressions of the classi¬ 
cal, neoclassical and turbulent particle fluxes are included 
in (Pq )®, (Pa^)®, and (P^respectively. The latter two 
fluxes are evaluated by the solutions (E'ai)ens and ha of 
the first-order drift kinetic and gyrokinetic equations, re¬ 
spectively. 


B. Ensemble-averaged energy balance equation 

The ensemble average of the energy density defined by 
Eq. (IMl) is written as 

{£)ens = I ^-0^-0 + ^+0{S), (93) 

a 

where the energy density of the electric held is neglected 
as a small quantity of O(d^). It is shown in Ref. 18 
that the radial components of the hrst two terms on the 
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right-hand side of Eq. (1651) are double-averaged over the same way as mentioned after Eq. (l90l) for the classical 
ensemble and the flux surface to give particle flux (E^)®. 


^ rfC// DaFa (^Hysc) _ X • Vs 


= E 




0(6^). 


Here, the radial particle flux is given by 

Eqs. (I87l) - (l89l) and the radial heat flux {qaY is written 
as 


U^^^r = ici^^r + i<iar ( 95 ) 


which consists of the nonturbulent part, 

( 9 ^)^ = (^JduJd^,J d^ Da{F,Yens{^^r^)e 

(I 5 

X f -niaU'^ + fJ-Bo - -TaO 


Vs 


= TaO [b X (v • ©a®"^)] • Vs^ + 0{SY. (96) 


, Ca-Sg L 

and the turbulence-driven part 


(9a )' = -{{^JdUjdfljd^ Daha{VYa X b) • Vs 


— iriaU'^ + IJ^Bq — 2^00 


0{SY- (97) 


In Eq. (IMl) . the heat stress tensor ©a^^ is defined 
by Tao©^°^ = JdUjdfiJd^DYFai)ens + 

fiBo — ^Tao)[maU'^hh + ^i?o(I — bb)]. The expression 
of Eq. (|M)) coincides with that of the neoclassical radial 
heat flux in terms of the heat stress tensor.— The tur¬ 
bulent heat flux in Eq. (EZl) takes the same form as that 
given by the conventional gyrokinetic theory.— 

The radial component of in Eq. (l69ll is ensemble- 
averaged to yield 


((QC . Vs)) = {q^r + ^rao(rC)« + o{6Y, (98) 


where the radial classical heat flux is given by 


( 9 ?r = E 


Sa-^O 


dU J d^ J F>a[{v X b) • Vs] 


X [CPMal,FtM)+C^,YFaMjbl)] 
^niaU'^ + fiBo - ^TaO 


= T 


oO 




■(Fa 2 X b) ■ Vs ) . 


(99) 


Here, Fa 2 = / d^v(maU^/2Ta — 5l2)ma'vCl is the col- 
lisional heat friction. The expression of the classical 
heat flux (q^)^ in Eq. (1991) agrees with the conventional 
one^ and it immediately gives the classical heat trans¬ 
port equation relating (q^)^ to the gradient forces in the 


Now, Eq. (155)) is rewritten as 


(94) 


at 


V' 


9 E ^‘‘oFaO ' ° 


Stt 


ds 


V' 


E + ^^-o(ra)y + • Vs) 

n X / 


- ( ^ E ^^oFaO + (Us • Vs) 


= dU J dfj. J d^ DaSa (^^rUaU'^ + hBq^'^ 

+ OiSY, (100) 


where the total radial heat flux is given by the sum of 
the classical, neoclassical, and turbulent parts as 


(9a)* = (9i®^))* + [q^y 

= iqa^r + {qtr + {qa)y (loi) 


and = (c/47r)(E)ens x Bq represents the non¬ 

turbulent part of the Poynting vector. Using the relation 
{d{Bl/iTT)/dt) = -(U')”^i9(U'(S(P°y"‘”s) • Vs))/ds - 
(Jo ■ (E)e„s) shown in Ref. 18, we also obtain 


'Wt 






E + ^Tao(ra)*] 

n. X / 


- ^ E ”“ 0 Ta 0 (Us - Vs) j 

= V'{3o-(E)ens)+V'Y,( f dU 

a ' 

X ^-niaU^ + AiHg^ ^ -(- 0{S^). 


dfi 


d^ DaSa 


( 102 ) 


Equations (I100|) and (11021) take the well-known forms of 
the energy balance equations^ except that the terms as¬ 
sociated with the electric field energy and the kinetic en¬ 
ergies due to the fluid velocities are neglected here as 
small quantities of higher order in 5. 
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C. Ensemble-averaged toroidal angular momentum 
balance equation 


The ensemble-averaged toroidal angular momentum 
balance equation is written as 


m 


V' 


'^naomaiuaiih + Ub) 


g(Poynting) 


■ec 




^{(nNA)^ + (nAr + (nC) 

a 

g(Poynting) 


riaomaiuaiih -I- Ue) + 


- — ((Vs 
47r 


e^(us • Vs) 

•eel 


ElEl 4- BB -I- (V X B)A 
= dU J d^,J D,Sam,Ub^^ + 0(5^), (103) 


where Ua\\ represents the nonturbulent part of the parallel 
fluid velocity for particle species a defined by naoWa|| = 
J dU J d^ij d^(Fai)ensU and = c(E)ens X h/Bo is the 
nonturbulent part of the E x B drift velocity. Equa¬ 
tion (I103D is derived from Eq. (1771) following the same 
procedures as shown in Ref. 18 except that the additional 
transport flux IlJ defined in Eq. (110711 and the external 
momentum source are newly included in the present case. 

On the left-hand side of Eq. (11031) . the terms including 
(Ma||b -4 Ue) and are of 0{6^) although they 

are written down to explicitly show the inertia-term part. 
The nonturbulent and turbulence-driven parts of the ra¬ 
dial flux of the toroidal angular momentum are defined 
by 


(n^A)^ = (^JdU JdflJ d^Da(Fal)ens 

X maC/6^ )ens ■ VS^ , (104) 

and 


(n^)^ 



dfi 


Daha 


X ma(Uh + Vao±) ■ • Vs)^^ , (105) 

respectively. It is shown in Appendix D that Eq. ((7^ is 
ensemble-averaged to give 


(((nC)«)) = ^nc + o(<53), (106) 

a 


where the radial transport flux of the toroidal angular 
momentum for species a due to the collision term and 
finite gyroradii is defined by 


= 


0 

^ [C'Qj((E'al)ens, EhM) + C^^(FaM, (E(,i)ens)]) ■ 


( 107 ) 


The expressions for the toroidal momentum fluxes shown 
in Eqs. (I104I) - (I107D agree with those given by conven¬ 
tional recursive formulations in Refs. 33-35. [Since the 
so-called high-flow ordering is used in Refs. 33 and 34, 
the expressions for the toroidal momentum fluxes in it 
reduce to those in the present work in the low-flow-speed 
limit.] As argued in Refs. 18 and 35, when there ex¬ 
ists the up-down symmetry of the background magnetic 
field, all toroidal momentum fluxes vanish to 0(6'^) and 
the nontrivial toroidal momentum balance equation is 
of 0{S^). In this case, gyrokinetic systems equations 
of higher-order accuracy in S are required for the cor¬ 
rect derivation of this 0{S^) toroidal momentum balance 
equation to determine the profile of the radial electric 
field^ although we should note, at the same time, that 
the radial electric field is not necessary to determine the 
particle and energy transport fluxes to the lowest order 
in 


VII. CONCLUSIONS 

In this paper, particle, energy, and toroidal momen¬ 
tum balance equations including collisional and turbu¬ 
lent transport fluxes are systematically derived from the 
gyrokinetic Boltzmann-Poisson-Ampere system of equa¬ 
tions. Considering an imaginary collisionless system, 
for which the distribution functions and electromagnetic 
fields coincide instantaneously with those for the consid¬ 
ered collisional system, and expressing the variation of 
the action integral for the collisionless system in terms of 
the solution to the governing equations for the collisional 
system clarify effects of the collision and external source 
terms on the collisionless conservation laws derived from 
Noether’s theorem. The gyrokinetic collision operator 
is newly presented, by which the collisional changes in 
the velocity-space integrals of the gyrocenter Hamilto¬ 
nian and the canonical toroidal angular momentum can 
be written in the conservative (or divergence) forms. It 
is confirmed that, to the lowest order in the normalized 
gyroradius, the ensemble-averaged fluxes in the derived 
particle, energy, and toroidal angular momentum balance 
equations can be written by the sum of conventional ex¬ 
pressions of classical, neoclassical, and turbulent trans¬ 
port fluxes. The extension of the present work to the 
case of the high-flow ordering remains as a future task. 
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Appendix A: COORDINATE TRANSFORMATION 

We consider the transformation of the phase-space co¬ 
ordinates in this Appendix, where the subscript repre¬ 
senting the particle species is omitted as far as it is 
unnecessary. In terms of the position x and the veloc¬ 
ity V of a given particle, we define the parallel velocity 
uy = v-b(x,t), the perpendicular velocity vj^ = v —riyb, 
and the zeroth-order magnetic moment, 


An arbitrary scalar field A on the phase space can be ex¬ 
pressed in terms of either the gyrocenter coordinates Z = 
(X, U,or the particle coordinates z = (x, r;||, po, ^o) 
as 

A^Z)=AP(z). (A5) 

Using Eqs. (IA4)) . (IMI) . and the Taylor series expansion, 
we obtain 


mv\ 

2Bo{x,t) 


(Al) 


where the equilibrium field at position x and time t is 
denoted by Bo(x, t) = i?o(x, t)b(x, t). We also define 
the zeroth-order gyrophase by = tan“^[(v • ei)/(v • 
82 ) where (ei, 62 , b) are unit vectors which form a right- 
handed orthogonal system at (x, t). Then, the gyrocenter 
coordinates Z = (X, t/, p,, are represented in terms of 
the particle coordinates z = (x, t!||,/io, 5 o) as 


AP{z) = (r*AP)iz) = AP{T{z)) = AP{z + Az) 


Es E A*"--A*'" 


n\ 

n—Q %\ 


5"A9(z) 


(A 6 ) 


where T* A^ denotes the pullback transformation of A® 
by T. Using the inverse transformation we also 

have A9(Z) = (J-^*AP){Z) = Ap{J-^(Z)). 

The Jacobians DP and for the two coordinate sys¬ 
tems z and Z are related to each other by 


(X, U, n, C) = (x, 1 ^ 11 , Mo, Co) + (Ax, Auji, A/ro, A^o), 

(A2) 

where 


Ax = —p -f O(S^) 

Az;|| = —?^||b • Vb • p — ^(3p • Vb • vj_ — vj_ • Vb • p) 
H-^i|| + 

rr> r> " 


A 


’0 


V , 

-^p ■ VBo 




+ 0{5^ 


= ^(v,.VB„)-gjp-Vb.p 


,_^vx-Vb-vx + jlL(b.Vb-vj 

n 


e 


g 

p ■ Vvi, ■ p H -• (Al X v) 


The formulas for A?;|| A^o, and A^o in Eq. (jA3l) are 
obtained by combining the guiding center and gyrocen¬ 
ter coordinate transformations Here, effects of the 

background electric field and turbulent electromagnetic 
fields are included through ij; [see Eq. m] and Al. 
When the background electric field and turbulent elec¬ 
tromagnetic fields vanish, Az;|| Apo, and ACo in Eq. (IA3I) 
agree with the results in Ref. 37. 

Denoting the coordinate transformation by T, 
Eq. (IA2I) is rewritten as 


(A4) 


DP{z) = det 


' djzy 

.d{z)_ 


D®(Z), 


(A7) 


where d{Z)/d{z) denotes the Jacobian matrix. Then, we 
use the following formula. 


z-Z) 
(AS) 

and partial integrals to derive the relation between the 
expressions of the scalar density BA in the gyrocenter 
and particle coordinate systems as 


^«(z+Az-Z) = f;l y: 


n\ 

n=0 ii , 


D®(Z)A®(Z) 

= J (fZ' J®(Z' - Z)D®(Z')A®(Z') 
= J (fz 5^{z + Az- Z)DP{z)AP{z) 


^ { -I Nn _ 

E 


n—Q 


5” [Az*i • • • Az^''DPlz)AP{z)\ 




z=Z 

(A9) 


where the replacement of z with Z is represented by 
[■ • -Jz^z = / (5®(z - Z) • ■ •. 


Appendix B: COLLISION OPERATOR IN GYROCENTER 
COORDINATES 

We can regard the collision term as a scalar field C on 
the phase space. When using the particle coordinates, we 
represent the collision term for collisions between species 
a and b by A well-established collision operator 

C^j(/a,/b) for the particle distribution functions fa and 
fb is known as the Landau operator [see, for example, 


Z = T(z) = z J- Az. 
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Eq. (3.22) in Ref. 32]. Then, the collision term rep¬ 
resented in the gyrocenter coordinates is related to 

by 


CliFa. F,) = r-^*ci,{r:Fa, TCFt), (Bl) 


and finite gyroradii of particles is defined by 

j5a(X) 

^ (-1)” V" 3" 

~ ^ (n-f 1)! ^ 

n—0 ^ ,• • • >'^Ti 


where the distribution function for species a {b) in the 
particle coordinates is written as the pullback fa = Fa Fa 
ifb = TffFb) of that in the gyrocenter coordinates Fa 
(Fb) by the coordinate transformation Ta {%) described 
in Appendix A, and 7^”^* transforms the collision term 
as a function of the particle coordinates into that of the 
gyrocenter coordinates. 

In order to see collisional effects on conservation laws, 
it is convenient to represent the collision term in the gy¬ 
rocenter coordinate using the transformation formula for 
the scalar density DaCab rather than that for the scalar 
Cab shown in Eq. (ED. Using Eq. (1X91) . we can derive 




k:=X 


d\ AxaC'P(z)AP(z) 


(B5) 


J x=X 


The integral of an arbitrary scalar field Aa over the whole 
phase space is written in either the gyrocenter or particle 
coordinate system as 


I FZDUZ)C^a{^)Al{Z) = I Fz Dl{z)Cl{z)Al{z). 

(B6) 

For the case of Ao = 1, Eqs. (IB4I) and (IB5|) reduce to 


Dl{Za)Clb[Fa,Fb]{Za)Al{Za) 

(-1)^ 


= E 


n=0 


n\ 


E 


9" [Azl^ ■ ■ • AzFDP{Za)C^ab[fa^ /fc] (z.) Ag (z, )] 
dzi^ ■ ■ ■ dzl^ 




(B2) 


where Aa is an arbitrary scalar field depending on parti¬ 
cle species and fa = T* Fa is rewritten by using Eq. (IMl) 
as 


faiZa) = 






n=0 


d'-^Faizg) 

dzl^---dz^ 


(B3) 


Then, the gyrocenter representation of the collision op¬ 
erator Cl^b acting on Fg and Fb is obtained by Eq. (IB2I) 
with putting A® = A^ = 1 and using Eq. (IB3I) to express 
fa and fb in terms of Fg and Fb, respectively. Integrating 
Eq. (IB2I) with respect to (U, fi, and taking the summa¬ 
tion over species b yield 


dU df, di DHZ)Cf,{Z)AHZ) 


d\ Cl{z)Al{z) 


- V-J 


Aa^ 


J z=Z 


(B4) 


where C® = Yl,b ^ab V = d/d'K. are used and 
J d^v = f dun / d^o J dfo DP{z) denotes the velocity- 
space integral using the particle coordinates. Here, the 
transport flux of the quantity Aa due to collisions 


I dU I d^b I d^ Dl{Z)CPa{Z) = -V • rC(X), (B7) 

and 



respectively, where f d^v C'^(z) = 0 is used. Here, is 
regarded as the classical particle flux which occurs due to 
collisions and finite gyroradii. In fact, using Axq ~ — p^, 
we see that the primary term of F^ shown in the last line 
of Eq. (jB^II is identical to the conventional definition of 
the classical particle flux F];* = (c/eaBo)Fai x b, where 
Fai = f d^vuiavCf is the collisional friction force. Thus, 
we have F^ = F^'[I + 0((5)]. 

Let us take the kinetic energy of the particle as Aa 
and put APg = ^mgvl = -I- poaHo(xQ). Then, it 

is written in terms of the gyrocenter coordinates as 

-^a = (^mavi^ 

= ^ITT-aUa + fiaBo{Xa) + HaPa ' VHo(Xa) 

‘^aUa{^V\\a'}za — 'Zta (^/^Oa)za—Z q, ^0 (^a ) H” ' ‘ ' 

= ^rUgUl + PaBoiXa) + ea{f^a{Za))u ” e„ (T"'* <^) (Z,) 

+ 0(^2), (B9) 

where the inverse 7])”^ of the transformation Fa given by 
Eq. (IA3I) is used. In this case, taking the summation 
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of Eq. (HH) over species a and using the conservation 
property Ea / rriaV^ = 0, we have 

^JduJdfiJ Di{Z)Cl{Z)T-^* 

= -V-Q^, (BIO) 

where represents the transport flux of the total ki¬ 
netic energy due to collisions and finite gyroradii defined 
by 


_ f -I '\n _ 


(n-kl)! ^ dXi^-'-dXir. 

a n=0 ^ H 


X J dV AxaAx^^ ■ ■ ■ Ax^j'CP(z)^maV^ 

= J2 J AxaCP{z)^ma'' 


J X=X 

(Bll) 


J x=X 


To the lowest order in 6, the collisional energy flux Qp is 
approximately written as Q® ~ Ea(‘ 3 a + Here, 

the classical heat flux for species a is defined by = 
{cTa/eaBo)Fa 2 X b, where Fa 2 = J d^'v{mav‘^/2Ta - 
5/2)toqvC'^ is the collisional heat friction. We note from 
Eq. (IB9I) that the expression of the kinetic energy in 
the gyrocenter coordinates should be generally given by 
the infinite series expansion in 6 in order for the gyro- 
center velocity-space integral of the collisional rate of 
change in the kinetic energy to take the form of the di¬ 
vergence of the energy flux without any local source or 
sink terms. In fact, this energy conservation property is 
broken if we keep only the lowest order terms in Eq. (IB9I1 
and evaluate the gyrocenter velocity-space integral 
J:JdUJd^^Jd^DS{Z)Ca,{Z){^m,U^ + ^,Bo{X)). 

The above-mentioned subtle relation between expres¬ 
sions of the collisional energy conservation properties in 
the particle and gyrocenter coordinate systems is also 
found when considering the collisional momentum con¬ 
servation. It should be recalled that the perturbative 
expansions in S are truncated up to finite orders in 
deriving gyrokinetic equations as shown in Sec. Ill al¬ 
though the conservative form of equations for the en¬ 
ergy and the toroidal angular momentum are obtained 
even from these approximate equations for the colli¬ 
sionless case since they are constructed based on the 
variational principle. Thus, from the viewpoint of 
practical applications, it is desirable for the approxi¬ 
mate collision operator in the gyrocenter coordinates 
to keep the conservation properties. More rigorously 
speaking, we want the gyrokinetic collisional velocity- 
space integrals Ea / dU J d^ J d^ D^{Z)C^{Z)Ha{Z) 
and EaJdUjdfiJd^ Ds{Z)Cs{Z){pl)s{Z) to take the 
divergence forms and include no local source or sink 
terms where Ha{Z) and (p^)®(Z) = (ea/c)H*^(Z) are 
the gyrocenter Hamiltonian and the canonical toroidal 
angular momentum defined by Eqs. (ES) and GH) , re¬ 
spectively. Here, it should be noted that not only kinetic 


parts of energy and toroidal momentum but also contri¬ 
butions from scalar and vector potentials are included in 
Ha{Z) and (ea/c)H*^(Z). In Appendix C, we find how to 
construct the approximate gyrokinetic collision operator, 
by which the two integrals mentioned above are written 
in the divergence forms. 

We now consider the entropy per unit volume de¬ 
fined in terms of the gyrocenter distribution func¬ 
tions as S'® = -J2al dU J dfj.J d(D3{Z) log Fa{Z), 
in which the rate of change is given by dS^/dt = 
-EafdUfdpfd^ Dl{Z) [logFjZ) + 1] {dFjdt). 
Then, the rate of change in S® due to collisions is ob¬ 
tained by putting A® = — [log Fa (Z) -|- 1] in Eq. (IB4II and 
taking the summation over species a as 


d^ Dl{Z)Cl{Z)[\ogF,{Z) + l\ 


= -E 


d^^Cliz) log/a(z) 


-V-Jg, (B12) 


<;=X 


where f d^v Cp(z) = 0 is used although we should recall 
that / dU f dfij d^F®(Z)C®(Z) does not vanish gener¬ 
ally as seen from Eq. (IBTl) . It is well-known that, when 
Landau’s collision operator is used for C^, the collisional 
entropy production rate given by the first term on the 
right-hand side of Eq. (IB12|) is nonnegative. This is 
Boltzmann’s H-theorem which proves the second law of 
thermodynamics. The collisional transport flux Jg of the 
entropy in Eq. (jB12l) is defined by 


J§(X) 

(n-fl)! ^ • • • 9A*" 

a n—0 ^ ' iiF" >'^n 


AXaA^i ■ • • Ax^" (7^(2) [log /a(z) I] 


x:=X 


= -E 


d^V AXaC'P(z)[l0g/a(z) 1] 


J x=X 


(BI3) 


It is shown that, to the lowest order in 5, the col¬ 
lisional entropy transport flux is written as Jg = 
Ea(‘^aoUa -f /To) where the lowest-order entropy 
density Sao for species a is given in terms of the lo¬ 
cal Maxwellian distribution function FaM as Sao = 
— f dU f dpf d^ FaM log FaM, and u))' is defined by 
Ua = Fa/ria- Here, we note again that the infinite series 
expansion in S as given by Eq. (IB2I) is used in deriving 
Eq. (IB12I) . When the expansion is truncated to finite 
order, the collisional entropy production term is repre¬ 
sented by — Ea I d^'v C® log fa plus residual error terms 
of higher order in d, and thus the H-theorem is only ap¬ 
proximately satisfied. 
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Appendix C: COLLISION OPERATOR RELEVANT FOR 
GYROKINETIC CONSERVATION LAWS 


In this Appendix, we consider an approximate gy- 
rokinetic collision operator instead of the one given by 
Eq. (IBII) [or Eq. (IB2I) with A® = = 1] in order to get 

the gyrokinetic collisional velocity-space integrals of en¬ 
ergy and canonical toroidal momentum to take desirable 
conservative (or divergence) forms. The approximate col¬ 
lision operator is written in the gyrocenter coordinates as 


OO - 


.n—0 

X E 

2 l ,••• , 2 , 

_ 

dXa 

d 


dvn 


dxa ■ ■ ■ dx]r 

[fa JbliZa) 

[Av^\aDPMC:dfaJb]M] 




[A/i 0 o 71 o(Za)C'^b[/o, /b](Zo)] 


J z.a=Zal 


(Cl) 


where Ar;||(j and A/roa are written as 

Axa = -h Ax(.2), Az;||„ = -F Az;j^), 

A^oa = A4'j+A42j, (C2) 

and fa is given from Fa by fa{za) = Fa{xa -f Ax^, ?;||oa -f 
Avlla,^J. 0 a + A/roo). Here, Aw[|^\ and A/i^^^ are the 0{S) 
parts of Aw||a and A/iga given in Eq. (IA3I) . In this Ap¬ 
pendix, we do not derive expressions for the 0{d^) parts 
Axa , AtiJi^, and AfigJ by the Lie perturbation ex¬ 
pansion method which is used to define the gyrocenter 
coordinates with the well-conserved magnetic moment 
because it would unnecessarily give higher-order accu¬ 
racy to the coordinate transformation than the accuracy 
of the gyrocenter motion equations themselves shown in 
Eqs. dSlli-dSIl). Instead, we determine these 0(6'^) terms 
from the conditions that the collisional change rates of 
energy and canonical toroidal angular momentum per 
unit volume in the gyrocenter space can be given in the 
conservative forms as shown below. Thus, the O(S^) 
terms are introduced not for accuracy of higher order 
in S but for satisfying the conservation property of the 
collision operator. 

In Eq. (ICID . the expansions in (Axi^\ Auya, A^Oa) are 
truncated to the first order while the infinite series ex¬ 
pansion in Axi^^ = —p^ is retained because fluctuations’ 
wavelengths in the directions perpendicular to the equi¬ 
librium magnetic field can be of order of the gyrora- 
dius Pa- In the WKB (or ballooning) representation, the 
above-mentioned infinite series expansion can be treated 


using the Bessel functions of the gyroradius normalized 
by the perpendicular wavelength i^i^i^° We should also 
note that the gyrophase average (• • • is taken so that 
the gyrokinetic equation with the collision term is solved 
only for the gyrophase-averaged part of the gyrocenter 
distribution function. 

For an arbitrary function A®(Za) which is independent 
of the gyrophase ^a, we obtain the following formula. 


/ r 1 

= ( DP{Za)CP,,[fa, MMA^M + E - 

\ L ”■ 

y- ■ ■ ■p-"i^g(za)Crj/a,A](z.)Ag°(zJ] 

"" dFf---dx:r 

- ■ ^^a^ DP{Za)CPf,[fa, fb]{Za)Al{Za) 

- [Au||aT>g(Za)C'^J/a,/b](Za)AS(Za)] 

- [A/^OaT>S(Za)C'PJ/a,/b](Za)A®(Za)] ) , 

A'L^^’Lal 

(C3) 

where A^iza) and AP°(za) are defined by 

^a°M = AliXa- Pa,Vlla,POa)- (C4) 

We should note that the function AP{za) defined from 
A®(Za) in Eq. (IC4I) does not exactly coincide with that 
given in Eq. (IA6I) in Appendix A by the second- and 
higher-order terms in the series expansion with respect 
to Azo. Integrating Eq. (IC3I) over the gyrocenter velocity 
space, we immediately obtain 


dU dp dei^®(Z)C®(Z)A®(Z) 


d^vCP{z)APM 


- V-J 


Aa'i 


J x=X 


(C5) 


where the transport flux due to collisions and finite 
gyroradii is defined by 

JAa(X) 


E 


{-Pa + Axi2)]C'P(z)AP(z 

-1 9 " 


J x=X 


(n + 1)! 


E 




p,p:f---p:-cp{z)AP{. 


■ (C6) 


J x=X 
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To the lowest order in S, Eqs. (IC3I) . (IC5I) . and (IC6I) de¬ 
rived from the approximate collision operator in Eq. (EH) 
agree with Eqs. (IB2I) . (IB4I) . and (IB5I1 given in Ap¬ 
pendix B, respectively. The particle flux due to col¬ 
lisions and finite gyroradii is given from Eq. (IC6I) with 
putting AP(z) = 1 in the same way as in Eqs. (IB8I1 . 

Now, let us take A®(Z) = Ha{Z) in Eq. (IC5I) . Here, 
Ha{Z) denotes the gyrocenter Hamiltonian defined by 
Eq. d^ni). It is desirable that the gyrocenter velocity- 
space integral J2a I dU J dfj. J D3{Z)CS{Z)HaiZ) 
takes the conservative form, which implies that the in¬ 
tegral is expressed by the divergence term only and 

^ f d^v Cl{z)Hl{z) = 0 (C7) 

a 

holds. Here, using Eqs. (E51) and (IC4I) . HP{z) is given by 


where / dU f dfj, f d^ D^{Z)CP{Z)Ha{Z) is not com¬ 
pletely given in the conservative form, the residual term 
Tha f d^'v CP{z)HP{z) = 0{d^) is smaller by a factor of S 
than other transport terms of C>(S^) in the lowest-order 
energy balance equation given by Eq. (llOOl) in Sec. VLB. 

We next put A®(Z) = (p^)®(Z) in Eq. (jC5p . Here, 
(p^)P(Z) denotes the canonical toroidal angular momen¬ 
tum defined by 

(p^)®(Z) = ^AJ(Z) = ^Aoc(X) + mMbdX), (Cll) 

where = — y and = I/Bq- We now see that 
EaJdUjdpJd^ i4f(Z)C'f(Z)(p-)9(Z) takes the con¬ 
servative form if 

Y, fd\CPdzWcrM=0. (C12) 




1 

2 


maUy„-|-/rHo(Xa)-|-ea(?!>(Xa)-|-AiJa(Za), (C8) 


where 

AHaiZa) = f Ax^2) . ^ + A^OaW^) Ha{Za) 

\ aXa dv\\a OfiQa J 

+ MOa [Bo{Xa - Pa) “ Bo(Xq)] 

-k Ba [^-IXa - Pa, V\\a, POa) “ (^(Xa)] . (C9) 

It is easily seen that Eq. (IC7|) is satisfied if AHa{za) = 
0. Then, substituting Eq. (IC2|) into Eq. (IC9I) and using 
AHa{za) = 0, we have 


Axl^) 


9xq 




Ik Qy 


Ik 




dpoa 


Ha{Za) 


— /^Oq 


So(Xa - Pa) - Bo{yia) + {Pa ' VHq) 1 -k 
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We find that the right-hand side of Eq. (ICIOI) is of 0{S‘^). 
Then, as remarked after Eq. (IC16I) . we can choose Axi^\ 
Au||^\ and Ap^J which satisfy Eq. (jC10|) and are of 
0{S^) so as to be consistent with Eq. (IA3I) . 

When we use Ax^^) = 0, Au||a = Au||]j\ and Apoa = 
foi' Eq. jCH) by putting Au[|^^ = Ap^^J = 0, we 
have AHa{za) = 0{5d and X)a/ d^^ C'P(z)i7P(z) = 
0((5^) because C'p(z) = 0(S) holds for the distribu¬ 
tion function, the zeroth order of which is given by 
the local Maxwellian. Therefore, even for this case 


Here, using Eqs. (IC4l) and (ICllI) . (p^)P(z) is given by 


(pVa(^) = ^[^0C(X) + Ai(;(x)] -k TOal^C + MPi)a(z), 

(C13) 

where 


A(p^)a(Za) = ^Ax^2)^ -k Au||a^^^ (x^ , Uya) 

+ (pdai^a - Pa,Vlla) 

- ^[Aoc(Xa) -k Ai^(Xa)] - rUaV^a- (C14) 

Again, we easily see that Eq. (IC12I) is satisfied if 
A(p^)a(za) = 0. The O(S^) variables, Axi^^ and Au|i^\ 
which meet the condition that A(pda(za) = 0, are given 
by 
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and 

7j|-bc)a(^a) = ^abdXa)Avj^l^ 

^^||a 

= -maVlla [bd^a “ Pa) “ bcM + Pa ' V^cl 
+ -f X(Xa - Pa) - X(Xa) + Pa ' - ^PaPa 


: VVx , 
(C16) 


where V = d/dxa and = —(VR-Vx)/{RBo) + ^bi^h- 
(V X b). Asa solution to Eq. (IC15(1 . we can assume Axi^^ 
to be given in the form Axi^^ = Aii^Vy. We should 
note that Vy • V(p^)® = 0{5~^) and 9(p^)®/9u||a = 
mab(^ = 0[5^) while the right-hand sides of Eqs. (IC15I) 
and (IC16I) are of 0{6) and 0{6‘^), respectively. There¬ 
fore, Eqs. (IC15I) and (IC16I) give Axi^^ and Au||^\ which 














































19 


are both of 0{S'^), consistently with Eq. (IA3I) . Then, 
these and Aw||^^ are substituted into Eq. (ICIOI) to 

determine of 0(6'^) as well. 

Thus, the collision operator, which has the desired 
conservation properties as well as the accuracy required 
for correct description of collisional transport of the en¬ 
ergy and the toroidal angular momentum, is given by 
Eq. (jCll) . in which Au||a, and A^oa are defined by 

Eqs. (IC2D . (ICIOI) . (IC15I) . and (IC16I) . Using this collision 
operator, putting .A^(z) = ^maV^^+iioBo{x.)+ea4>{^) and 
AP{z) = (ea/c)[Aoc(x)-|-Ai^(x)]-|-maU^ in Eq. (IC6| and 
taking their summation over species a define the trans¬ 
port fluxes Q® and of the energy and the canoni¬ 
cal toroidal angular momentum, respectively, which ap¬ 
pear in the energy and toroidal angular momentum bal¬ 
ance equations in Secs. V.A and B [see Eqs. (IB51) . (IMl) . 
(1771), and (1751) ]. In the definition of mentioned above, 
the contribution of the potential energy part Cqc/) is writ¬ 
ten as which is smaller than the contribution 

of the kinetic energy part by a factor of S because the 
classical particle fluxes represented by the lowest-order 
part of are intrinsically ambipolar. Therefore, the 
energy flux Qp defined here agrees with Eq. (IBllI) to 
the lowest order in S. Regarding the entropy produc¬ 
tion discussed in Appendix B, the positive definiteness 
of the entropy production rate [corresponding to the first 
term on the left-hand side of Eq. (IB12I) ] is only approxi¬ 
mately shown by using the present model collision oper¬ 
ator in Eq. (IC5I) with .4.®(Z) = — [logFQ(Z) -|- 1] because 
AP{z) = — [log/a(z) -I- 1] is not rigorously derived from 
Eq. (IC4I) without the infinite series expansion in Az^ as 
in Eq. (lAB . 


Appendix D: Derivation of Eqs. (I106II and II107II 


In this Appendix, it is shown how to derive Eqs. (11061) 
and (11071) by using the collision operator given in Ap¬ 
pendix C. On the right-hand side of Eq. (1751) where the 
radial flux (B®)® of the toroidal angular momentum due 
to collisions and finite gyroradii is defined, the two types 
of fluxes and evaluated 

by taking the summation of Eq. (|C6I) over species a with 
putting AP{z) = Ca and Ap{z) = (p^)P(z), respectively. 
Here, (p^)S(z) = (ea/c)[Aoc(x) -h Ai^(x)] + maV(; is used 
for the collision operator which conserves the toroidal 
angular momentum as explained in Appendix C. Conse¬ 
quently, the ensemble average of (H^)® is expressed ex¬ 


plicitly up to 0{6^) as 

((nC)®)ens = ((j;rC + fji) 

= J2[ rf'v(C'P(z)) 
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(Dl) 


from which Eqs. (I106|) and (I107|) are immediately ob¬ 
tained. It is noted that the O(S^) part of ((n^)®)ens has 
no contribution from the gyrophase-dependent part of 
the distribution function, the lowest-order part of which 
is given by fai = —Pa ' '^FaM with the gradient operator 
V taken for the fixed energy variable e = imaU^-|-e((/))ens- 
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